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Abstract 



for the wave covectors. For the class of general nonlinear Lagrangian models, 



. We study the wave propagation in nonlinear electrodynamical models. Partic- 

00 . 

CN . ular attention is paid to the derivation and the analysis of the Fresnel equation 

o 

we demonstrate how the originally quartic Fresnel equation factorizes, yield- 

~o 

Q^. ing the generic birefringence effect. We show that the closure of the effective 

bJQi constitutive (or jump) tensor is necessary and sufficient for the absence of 



birefringence, i.e., for the existence of a unique light cone structure. As an- 
other application of the Fresnel approach, we analyze the light propagation 
in a moving isotropic nonlinear medium. The corresponding effective consti- 
tutive tensor contains non-trivial skewon and axion pieces. For nonmagnetic 
matter, we find that birefringence is induced by the nonlinearity, and derive 
the corresponding optical metrics. 
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1. INTRODUCTION 



Wave phenomena belong to the most interesting and important processes in physics. 
Among other field theories, nonlinear electrodynamics attracts much attention in connection 
with the prominent role played by light in the experimental and theoretical studies of the 
structure of spacetime and matter. 

Nonlinearities in electrodynamical models can arise in different ways in classical and 
quantum field theories. For example, the old Born-Infeld theory [ED] was a fundamental 
theory alternative to the classical Maxwell electrodynamics which provided a model of a 
classical electron. On the other hand, quantum Maxwell electrodynamics predicts nonlinear 
effects which arise due to the radiative corrections, see p|-§J. Finally, in the modern string 
theories a generalized Born-Infeld action naturally arises as the leading part of the effective 
string action, see for example. 

Wave propagation in the various nonlinear electrodynamical theories was studied previ- 
ously in f8HT2|| and also in fl3HT6| . A general feature revealed in these studies is the existence 
of birefringence. In crystal optics, the notion of birefringence means the emergence of two 
rays (ordinary and extraordinary) with different velocities inside the material medium. We 
will use the expression "birefringence" in a similar sense, associating it with the situation 
when two different light cones exist for the wave normal covectors. However, the earlier 
results are incomplete in the sense that the full Fresnel equation, governing the wave nor- 
mals, was never derived explicitly. Moreover, it was not demonstrated how it happens that 
the original quartic surface of wave normals reduces to the light cone. That is the primary 
interest in our study, and we will try to clarify this aspect for the nonlinear electrodynami- 
cal models, using and expanding our earlier results obtained within the framework of linear 
electrodynamics pl| , |23[| . 

Our basic tool will be the general formula for the Fresnel equation derived earlier within 
linear electrodynamics. Now we observe that the analysis of the wave propagation in a 
general nonlinear model reduces to the linear case because the jumps of the derivatives of 
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the excitation and of the field strength are in all cases related by a linear law. We can make 
then use of our master formula for the Fresnel tensor and derive the Fresnel equation for 
any nonlinear model, and thereby explain the reduction to the light cones. 



2. ELECTROMAGNETIC WAVES AND FRESNEL TENSOR 

Quite generally, Maxwell's equations for the excitation 2-form H = (T>, 7i) and the field 
strength 2-form F = (E, B) read 

dH = J, dF = 0. (2.1) 

Here J is the electric current 3-form. These equations must be supplemented by a con- 
stitutive law H = H(F). The latter relation contains the crucial information about the 
underlying physical continuum (i.e., spacetime and/or material medium). Mathematically, 
this constitutive law arises either from a suitable phenomenological theory of a medium or 
from the electromagnetic field Lagrangian. It can be a nonlinear or even nonlocal relation 
between the electromagnetic excitation and the field strength. 

If local coordinates x l are given, with i,j, ... = 0, 1, 2, 3, we can decompose the excitation 
and field strength 2-forms into their components according to 

H = - Hij dx { A dx j , F = - F i3 dx { A dx j . (2.2) 

We will study the propagation of a discontinuity of the electromagnetic field following 
the lines of Ref. see also Refs. |[23|| . The surface of discontinuity S is defined locally by 
a function $ such that $ = const on S. Across S, the geometric Hadamard conditions are 
satisfied: 

[i^] = 0, [d i F jk ] = q i f j k, (2-3) 
[H^ = 0, [diH jk ] = qi h jk . (2.4) 

Here [J-\ (x) denotes the discontinuity of a function T across S, and := <9j$ is the wave 
covector. Given the constitutive law H(F), which determines the excitation in terms of the 



field strength, the corresponding tensors and htj, describing the jumps of the derivatives 
of field strength and excitation, are related by f24 



1 dH 

hij = - Kij kl fki, with Kij kl :=— ^. (2.5) 

We will call k^ 1 the jump tensor. In linear electrodynamics, its components coincide with 



the components of the constitutive tensor (which describes the linear law = | F) 



kl), 

and they are independent of the electromagnetic field. However, in general the jump tensor 
function of the electromagnetic field, the velocity of matter, the temperature, and 
other physical and geometrical variables. Quite remarkably, all the earlier results obtained 
for linear electrodynamics remain also valid in the general case because whatever local 
relation H(F) may exist, the relation between the jumps of the field derivatives, according 
to ( p.5|) , is always linear. 

If we use Maxwell's equations ( |2.1|) , then (|2.3| ) and ( |2.4|) yield 

e ijkl qj h kl = 0, e ijkl qj f kl = 0. (2.6) 

Let us introduce the analog of the conventional constitutive matrix 

ijkl ± ijmn kl _ Ul (r) 7 \ 

where we denote 7i u := | e l i mn H mn . Similarly to Kij kl , we will often call the tensor x^ kl the 
jump tensor density. 

Now, making use of ( |2.5|) and Q2.7|) , we rewrite the system ( |2.6| ) as 

X ijkl qjfki = 0, e ^ kl qj f kl = 0. (2.8) 

Solving the last equation by = q^aj — qjcii, we finally reduce (|2.8|)i to 

X^qma^O. (2.9) 

This algebraic system has a nontrivial solution for only when the wave covectors satisfy 
a certain condition. The latter gives rise to our covariant Fresnel equation |pl] , |23 l 



G m (x)mq k qi = o, (2.10) 

with the fourth order Fresnel tensor density Q of weight +1 defined by 

G ijkl (x) ■= ^ mnpq i rstuX mnr{l x jlpslk x l)qtu - (2.11) 

It is totally symmetric, Q l ^ kl {x) = G^ kl \x)i an d thus has 35 independent components. 

3. NONLINEAR ELECTRODYNAMICS 

Let us denote the two independent electromagnetic invariants as 

/. : l-',,l-" J . h : I',, l" J . (3.1) 

where i™ = |^ 4jfci ' Fm and rf^ kl := \g\~ l / 2 e l ^ kl . The Hodge operator for the exterior forms is 
denoted by the star *, as usual; but we will use a tilde to denote the dual 2-tensors. We will 
not restrict ourselves to the case of Minkowski spacetime with = diag(l, —1, —1, —1), but 
instead will be considered as an arbitrary curved Lorentzian spacetime metric. 

The class of nonlinear electrodynamics models we study are described, in general, by the 
Lagrangian 4-form 

V = Lr), with L = L(I U I 2 ). (3.2) 

Here, as usual, rj is the 4-form of the spacetime volume. The electromagnetic excitation 2- 
form, which enters the Maxwell equation ( |2.1| ), is derived as the derivative of the Lagrangian 
form, H = —dV/dF. Explicitly, we then have the nonlinear constitutive law 

H = 4 (— L 1 *F + L 2 F) . (3.3) 

We denote the partial derivatives of the Lagrangian function w.r.t. its arguments as 

dL , d 2 L 



~ a ' dij ^ ab - di a di b ' 

In accordance with (|2.5| ) and ( |2.7|) , the direct differentiation of ( |3.3|) yields the jump tensor 



a, b =1,2. (3.4) 



k x g i[k g l]j + k 2 F 13 F kl + k 3 F ij F kl + h F ij F kl + k 5 F ij F kl + k 6 rf jkl \ . (3.5) 
The coefficients kA, A = 1, . . . , 6, are functions of the electromagnetic fields: 

k x = 4Li, k 2 = 8L U , k 3 = k A = 8L 12l k 5 = 8L 22 , K = 2L 2 . (3.6) 



The identifications ( |3.6| ) are derived for the nonlinear Lagrangian (|3.2|) from the constitutive 
law ( |3.3| ). However, in most computations below we will consider the most general case with 
unspecified arbitrary coefficients k^- This may be useful if we want to study the nonlinear 
electrodynamics of a more general type, for instance, with the dissipation effects and/or in 
moving media. 

In general, the untwisted tensor density x^ kl ( x ) °f weight +1 has 36 independent com- 



ponents. We can decompose it into irreducible pieces |23| with respect to the 6-dimensional 
( "bivector" ) linear group as follows: 

^ijkl _ (l)^,ijkl _|_ (2)^,ijkl _|_ (3)^ijkl ^ (5-7) 

The irreducible pieces of % are defined by 

(2) x ijkl ._ 1 ^ijkl _ ^klij^j = _ (2)^«ij ^ (3) x %jkl ._ x [ijkl] ^ ^.$) 
Wyijkl ._ ^ijkl _ (2)^ijkl _ (S)^ijkl _ m^klij _ g-j 



The irreducible pieces ^x, and ™x have 20, 15, and 1 independent components, 

respectively. The possible presence of an axion piece ^% was first studied by Ni p5 



whereas a constitutive law with an isotropic skewon was discussed by Nieves and Pal 
2|. 

In the Lagrangian models (|3.2j ), the effective constitutive tensor is automatically sym- 



metric, i.e. ™x — 0, which follows from (|3.6|) , since &3 = ^4. However, in general the jump 
tensor ( |3.5|) has all the three irreducible pieces: 



(l)yij'W 



fcl /yii + fc 2 F fc/ + ih±M. (p** F kl + F ij F kl 
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+ k 5 F ij F kl - -- [{k 3 + h)h + {h - k 2 )I 2 ] V 



ijkl 



12 



(2) ijkl 



(3) ijkl 



^ pij pkl _ pij pkl 



-^J\g\ Ps + h)h + (k 5 - k 2 )i 2 + i2fc 6 ] t^' w . 



(3.10) 

(3.11) 
(3.12) 



4. FRESNEL EQUATION AND BIREFRINGENCE 

Our study of the algebraic system in the framework of the Hadamard formalism yields the 
Fresnel equation in the generally covariant form (|2.10|) with the Fresnel tensor density ( p.ll|) . 



For the explicit jump tensor density (|3~5|) , it thus remains to substitute its components into 



( p. 11Q - A straightforward calculation yields the result: 



gijkl = _ g^g k ^ + iy g {ij t kl) + Zt {ij t kl) ) . (4.1) 



Here we denote 



f 3 := F F J k , (4.2) 



and 



X = k\ + ^ (k 3 + k 4 ) I 2 - hh h + \ (hh - k 2 k 5 ) if, (4.3) 

y = kt (k 2 + k 5 ) + {k 3 k 4 - k 2 k 5 ) h, (4.4) 

Z = 4 (k 2 k 5 - k 3 h) . (4.5) 



The most remarkable property of ( |4.1| ) is that it is obviously factorizable into a product 
of 2 second order tensors. Correspondingly, the quartic Fresnel surface of the wave normals 
reduces to the product of two second order surfaces: 

— ki —ki 

G l]H (x) QiWkqi = g^r (g\ J m) {g 2 q k qi) = {fl m) (gf qm)- (4.6) 

In other words, the wave normals lie not on the quartic surface but on one of the two cones 
which are determined by the pair of optical metric tensors: 



g? := Xg* + (y + Vy 2 - XZ) t», (4.7) 

gl j := X 9 v + Q? - Vy 2 - XZ) tK (4.8) 

The second equality in (4.6) offers a different description of the cones by means of the 
conformally equivalent metric tensors: 

H := (y - Vy 2 - XZ)g« + = ^(y~ Vy 2 - XZ) g{ j , (4.9) 



H : = {y + vy 2 - xz) g ij + zt ij = ^(y + vy 2 - xz) g <$. (4.10) 



Thus, the general Fresnel analysis demonstrates that in any nonlinear electrodynamics model 
( p.2[ ) the quartic wave surface always reduces to two light cones. This is the birefringence 
effect which is thus a general feature of the nonlinear electrodynamics. 



5. PROPERTIES OF OPTICAL METRICS 

Let us discuss the results obtained in the previous section. The following general obser- 
vations are in order. 

The Fresnel equation is trivially satisfied for all wave covectors when fci = 0, see fl4.ip . 
Thus, - in order to have waves - every electrodynamical Lagrangian L should necessarily 
depend on the invariant I\ = F^F 1 ^ (thus providing k\ ^ 0). 

Accordingly, we will always assume that fcj ^ 0. 

In order to have a decent light propagation, the optical metrics should be real and with 
Lorentzian signature. How can one be a priori sure that for every L an optical metric 
necessarily has these properties? 

Using ( |4.3|) -( f4~5f) we find an explicit expression for the quantity under the square root in 
the above formulas : 

y 2 -XZ = N 2 + N 2 N 3 , (5.1) 

where we have denoted 
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Ni := h (k 2 - h) + (hK - k 2 k 5 ) h, 
N 2 := 2hk 3 + (k 3 k 4 - k 2 k 5 ) I 2 , 
N 3 := 2hh + (hh - k 2 k 5 ) I 2 . 



(5.2) 
(5.3) 
(5.4) 



The expression ( |5.1[ ) is always non-negative in every nonlinear theory ( p.2| ) because N 2 = N 3 
when we take into account that k% = k^, see ( 13.61) . 

The signature of a four-dimensional metric is Lorentzian if and only if its determinant is 
negative. Straightforward computation yields for the optical metrics 



(det^) = (det^) 



a" + 



a(3 a 



h + sign(#) 



16 



a = 1,2. 



(5.5) 



Here a — X and Pi — y + Vy 2 - XZ, (3 2 = y - Vy 2 - XZ. As we see, both optical 
metrics have Lorentzian signature as soon as the spacetime metric is Lorentzian. 

Summarizing, we have demonstrated that fl4.7| )-( fL8|) indeed describe the generic effect 
of a birefringent light propagation for all nonlinear Lagrangians. 

Recently, the emergence of the two "effective geometries" has been described in [p~3|.[T~I 



without using the Fresnel approach. This result is in a qualitative agreement with our 
analysis. In order to prove the quantitative correspondence, one needs to show that our 
optical metrics are conformally equivalent to the effective metrics of []l3lJl4l. Although the 
corresponding comparison is rather complicated and the direct proof is still not available, 
one can verify that \/y 2 — XZ is equal to 64a/A of Jl^j. Moreover, one can recast the 
optical metrics (|4.7| )- ([4.10|) into the form of the so-called Boillat metrics of [ |I~5|j . 



6. SPECIAL LAGRANGIANS 

It is worthwhile to study in a greater detail certain particular nonlinear models which 
are potentially of physical interest. 
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A. Lagrangian L = L(I 2 ) 



When the Lagrangian depends only on the second electromagnetic invariant, we have 
k\ = 0, and there are no waves in such models. 

B. Lagrangian L = L(Ii) 

For the Lagrangian which, on the contrary, depends on the first invariant only, we find 
k 3 = k 4 = k 5 = 0. Accordingly, from (|4.3|) -(|4~5|) we find X = k\, y = k x k 2 , Z = and thus 
0-(H yield 

<g = h (h + 2k 2 t ij ) , $ = k\ g*. (6.1) 

Correspondingly, we still have birefringence with some photons moving along the standard 
null rays of the spacetime metric g lJ , whereas other photons choosing the rays null with 
respect to the optical metric Li g^ + 4Ln t 13 ' , cf. |T3|,p^ 



C. Lagrangian L = U{I\) + al 2 

This is a simple generalization of the above case. Here a does not depend on the electro- 
magnetic field, although it is not a constant, in general. When it depends on the spacetime 
coordinates, a = a(x), one can identify it with the axion field. 

Here we again have k^ = k 4 = A; 5 = and we recover the same light cone structure flS.l| ) . 
To put it differently, axion does not disturb the light cones which are solely determined by 
the spacetime metric and by the dependence of the Lagrangian on the invariant I\. 

D. Lagrangian L = al\ + V{I 2 ) 

Then k 2 = k 3 = fc 4 = which yields X — k\ — k\k§ I\, y = kik 5 , Z = 0. Consequently, 
gl j = h [(h - k 5 h) J* + 2k 5 t ij ] , g$ = h{h - h h) g ij , (6.2) 
i.e., there is again birefringence with one cone determined by the standard spacetime metric. 
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E. Born-Infeld theory 



The Lagrangian of the Born-Infeld (BI) theory |IJ reads: 



1 



2b 2 1 166 4 



II -1 



Here b is the coupling constant. By differentiation, we find: 
fci = 4Zq = ' 

fc 2 = 8Ln = 



1 + 1 r - 1 t^ 1 ' 2 

-1 



2^ 2 (l + ^A-T^/l) 372 ' 
k% = = 8-L12 



2b 2 16b 4 J 2 

- (1 + 2^1 



•22 



3/2 ' 



Correspondingly, 



X = 



y = 



z = 



(l + 2I2 A 16 & 4 



166 4 
1 

26 2 



b 2 {^ + 26 2 A leV I 2 



2 ' 



1 



& 4 ( 1 + M-1^I 

As a result, we find 3^ 2 — = 0, and thus the two optical 



(6.3) 



(6.4) 
(6.5) 
(6.6) 

(6.7) 



(6.8) 

(6.9) 
(6.10) 

metrics ( |4.7| )-(fO|) coincide, 



ij ij 
9l = 92 = 



1 I It 1_ 12 

1 T 2h 2 J l 1664 J : 



(6.11) 

26 2 11 166 4 J 2 

Birefringence disappears, and the photons propagate along a single light cone determined 
by the "quasi-metric" of Plebanski [jO]] 
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7. NO BIREFRINGENCE CONDITION 



In this section, we will restrict our attention to the Lagrangian theories for which the 
constitutive tensor is ( |3.5|) , and the coefficients are derived as (|3.6|) . It is important that the 
Fresnel analysis reveals that ki 7^ 0, otherwise there is no decent wave propagation at all. 

As it is clear from ( |5.1| ), the necessary and sufficient condition of the absence of birefrin- 
gence is provided by the pair of equations: 

Ni = h (k 2 - k 5 ) + (k 3 h - hk 5 ) h = 0, (7.1) 
N 2 = N 3 = 2hh + (hh - k 2 k 5 ) h = 0. (7.2) 

Here the property k% = ki of the Lagrangian models is used. 

Taking into account that all fc's are the partial derivatives of the Lagrangian L w.r.t. l\ 
and/or I 2 as displayed in ( |3.6| ), we can view the above system as a pair of partial differential 
equations, the solution L = L(/ 1 ,/ 2 ) of which describes a model without birefringence (i.e., 
with a single light cone). At least two such particular solutions are already known: one 
is rather simple, namely, the standard Maxwell theory with L = ii/4. Another is more 
nontrivial - this is the Born-Infeld theory with the Lagrangian (|6.3|) . One may ask the 



question: Are these the only solutions of the system (|7.1|) - (|7^2]) ? The immediate inspection 
of the system ( |7.1|) -( fT2|) shows that the answer is negative. For example, the Lagrangian 
function L(/ 1 ,J 2 ) = alxjl 2y with constant a, satisfies the equations ([71|)-( f7"2"D . Such a 
nonlinear (and nonpolynomial) model thus also has no birefringence. It is an open problem 
to find the complete set of solutions of (|7.ip-(f7"2"D, leading then to a single light cone. 

8. CLOSURE CONDITION 

Let us denote the "traceless" part of the jump tensor as 

^ijkl ._ (l)^ijkl _|_ (2)^ijkl^ (S.l) 

As we know |]23|| , only the traceless part determines the Fresnel surface, whereas the axion 
part ( 3 \ l i kl drops out completely from the wave propagation analysis. 
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For the general jump tensor (|3.5|) , we find: 

1 . 



4 e pqij 



it ^ ^klrs it 



+ 0,1 F pq F mn + a 2 F pq F mn + a 3 F pq F mn + a 4 F pq F 



Here the coefficients read: 



k\ 1 
oq = — [(h + h) h + (h ~ k 2 ) h] = 77 




a 1 



- (k 3 + fc 4 ) ki - fc 3 /c 4 1 2 + -j [2(k 3 + h) h + (2h + k 2 ) I 2 ] 



(8.2) 



(8.3) 



(8.4) 



(8.5) 



(8.6) 



a 3 = (k 5 - k 2 ) h + k 2 h h + [2{h - k 2 ) h + (2k 3 - fc 4 ) Ji] 

V 3h 7 2 ^ 7 ' 

a 4 = - (A; 5 - fc 2 ) fci - k 2 h h - — [2(h ~ k 2 ) h + (2h - h) h] 

The second lines in (|8.4j)-(|8~7|) give the a's in terms of the combinations ( |7.1|) -( [T2"D . Certainly, 
we use the assumption that ki 7^ 0. 



Like the constitutive tensor of linear electrodynamics, the jump tensor /t, 



ki 1 



2 ^-ijmn 



^rnnki determines a linear map in the 6-dimensional space of 2-forms. When the action of 
this map, repeated twice, brings us (up to a factor) back to the identity map, we speak of 
the closure property of jhij kl . The importance of the closure property is related to the fact 
that ultimately fa 1 * 1 turns out to the duality operator which determines a unique conformal 
Lorentzian metric on the spacetime. 

In nonlinear electrodynamics, the jump tensor ( |3.5|) has the closure property when 



a = 0, 
as is evident from (18.21). 



a\ = 0, ct2 = 0, 



a 3 = 0, a 4 = 0, 



(8.8) 
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9. EQUIVALENCE OF CLOSURE AND NO BIREFRINGENCE CONDITIONS 



In linear electrodynamics, there is much evidence (although the final rigorous proof is 
still missing) that the quartic Fresnel surface of wave covectors reduces to a unique light 
cone if and only if the constitutive tensor has the closure property. 

For the nonlinear Lagrangian theories it is possible to make some progress in solving 
the equivalence problem. In a certain sense, the situation here is simpler because we have 
discovered that the quartic Fresnel surface is always reduced to the product of the light cones 
(birefringence). The next step is thus to study under which conditions the birefringence 
disappears and, correspondingly, a unique light cone arises. 

Theorem: In the general nonlinear electrodynamical model described by the Lagrangian 
Q3.2|) , the Fresnel equation implies a single light cone (no birefringence) if and only if the 
traceless part of the jump tensor satisfies the closure property. 

Proof: As a preliminary remark, we note that for the Lagrangian models (|3.2| ) the jump 
tensor ( [3.5|) is symmetric because fc 3 = fc 4 . As a result, A 2 = A 3 . 

The necessary condition is evident. The birefringence is absent when Ai = N 2 (= A3) = 
0, see ( |7.2| ). Then we immediately read from (|8.3p -( ]57FD that a = a\ = a 2 = a 3 = a 4 = 0, 
and thus the closure is recovered from ( |8.2| ). 

The sufficient condition is also proved straightforwardly. The closure condition ( |8.8| ) has 
the unique solution 

N 2 + N 3 , s 

Ni = 0, — =— — - = 0, (9.1) 



when we analyze (|8.3| )-( p7f[ ). Since for the Lagrangian models we have A 2 = A3, then ( |9.1|) 
yields Ai = N 2 = A3 = 0. Thus, there is no birefringence. 

To put it differently, we have proven that the closure of the traceless jump tensor is the 
necessary and sufficient condition for the reduction of the fourth order Fresnel wave surface 
to a single light cone. This is true for all nonlinear Lagrangian electrodynamical theories 
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( p.2|) . Returning to our studies of the general linear electrodynamics, we expect that a 
similar result holds true there. 

A. On asymmetric jump tensors 

Symmetry of the jump tensor is very important in the equivalence proof above. In order 
to clarify this point, let us consider an arbitrary jump tensor Q3.5|) without assuming the 
explicit form of the coefficients (|3.6|) . When k 3 ^ k^, the jump tensor has the nontrivial 
skewon part Q3.ll]) . Also, N 2 7^ N 3 [in fact, as we can see from (|5.3| ) and Q5.4|) , N 2 — A 3 = 

We can easily verify that the closure of an asymmetric operator is not equivalent to the 
no-birefringence property. Indeed, take, for instance, k\ 7^ 0, k 2 = k 3 = 0, A; 4 7^ and 
&5 = 0. Then the jump tensor is asymmetric and Ni = N 2 = but A3 7^ 0. We then obtain 
a unique light cone because ( |5.1| ) vanishes identically. However, the closure condition ( |B.8| ) 
is not satisfied, since (|8.3|) - (j577D are nontrivial for A x = N 2 = and A 3 7^ 0. This also 
means that the requirement of a unique light cone does not necessarily implies that x must 
be symmetric. 

The opposite is also true: Suppose an asymmetric jump tensor ( |3.5|) has the closure 
property, i.e. (fg.8|) is fulfilled. Then we find ( ^Tl) again. However, (|9.1| b) yields A 2 = — A 3 , 
and consequently together with (|4.7| ) and (|4.8|) describe the case of a birefringent and 
dissipative wave propagation. 

These examples show that the closure of an asymmetric jump (or constitutive) tensor 
does not guarantee the absence of birefringence, and, vice versa, no-birefringence is not 
accompanied by the closure property for an asymmetric operator. 

10. MOVING ISOTROPIC NONLINEAR MEDIA 

Recently, there has been some interest in the light propagating in moving media with 
nontrivial dielectric and magnetic properties. The first covariant analysis of the Fresnel 
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equation for this case was done by Kremer [19). An isotropic medium is characterized by 
the constitutive law 



-F 11 +2[--e) u k F k[i u j] 



(10.1) 



where u l is the 4- velocity of the moving matter (normalized as usual by UiU 1 = 1), and e and 
H are the permeability and permittivity functions of the isotropic medium. The case when 
they do not depend on the electromagnetic field strength (being constant in space and time, 



for example) was investigated in [flQf| . 

More recently, the nonlinear case when e = s(F) and \x = f^{F) are functions of the 
electromagnetic field has been studied by De Lorenci et al [[L4|]. However, the attention was 
restricted to certain special cases, and the general result is still missing. 

We can perform a fairly complete analysis of the wave propagation in a nonlinear moving 



media on the basis of our covariant Fresnel equation ( 2.1C ). By differentiation, we easily 
find the jump tensor fl2.7|): 



X 



ijkl 



+ m kl F i3 +2(m kl - e kl ) u m F m[l u Jl 



(10.2) 



The second line is absent in the linear theory, and the tensors 



m 



dF, 



OF, 



10.3) 



are responsible for the nonlinear electrodynamical effects. 

Inspection immediately reveals that the jump tensor density ( |10.2| ) contains both an axion 
and a skewon part. There are claims in the literature that axion and skewon, in general, 
do not have physical sense. However, here we encounter a simple and a physically sound 
example with the both parts being nontrivial. The irreducible pieces (|3.8|) - (|3.9| ) read: 



jV [ y ]j + 4 Q - u [k g l][i u j] + ~ (m ij F kl + m kl F ij 
+ (rji _ e ^') PK 1 ^ + ^ rp< [F mn m mn + 2{m mn - e mn )P m u n ] 



10.4) 
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(2) ijkl 



- (F ll m kl - F kl m 11 ) + (m kl - e kl ) P [i u> ] - (m 11 - e ij ) P [k u l] 



10.5) 



and for the axion piece, we have 



ijkl 



\9\ e ijkl 



F ij m kl + 2 (m kl - e lj ) P 



V 



;io.6) 



Thus, nonlinear isotropic matter does have axion and skewon induced by nonlinearity. 



A. Nonmagnetic matter 



Let us consider the case when the magnetic constant is independent of the electromag- 
netic field, that is m v = 0. [In the simplest case, we can restrict the attention to the 
purely dielectric medium with \x = 1. However, we will formally keep /i ^ 1, for the sake of 
generality] . 

It is convenient to make the evident split of ( |10.2| ) into the sum y?i kl = </> y + , with 



<f) ijkl 



kl , 



2J\g\e kl u m F m[l u^. 



10.7) 
10.8) 



It is straightforward to find the Fresnel tensor for the first piece: 



Q l3k \<P) = sign(s) J\g] g& g kl \ 
H 2 v 



Here we have denoted the so called Gordon optical metric ]TB 



as 



10.9) 



gv ■= g v + (en - 1) U V 



(10.10) 



Its inverse reads 



( 1 A 

9ij ■= 9ij + I — -II UiUj, 

and the determinant can be easily computed: detg = (det g)/(e/j l ). 

Using this optical metric we can simplify the form of ( |10.7| ), bringing it to 



(10.11) 
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wki = * /|^| gi[kfii]j = 2,1- \/\g\ g i[k g l]j . 



(10.12) 



Now, for x = + -0, using a compact notation by omitting the indices, we have 



Q{x) = £(0) + W) + ^ (Oi + o 2 + o 3 + 21 + r 2 + t 3 ) . 



10.13) 



Here the mixed terms O a contain one ^-factor and the T a 's two ^-factors. Postponing the 
symmetrization over k, I to the very last moment, these terms read explicitly as follows: 



Ol(0, V,0) 


^mnpq ^rstu ' , 


(10.14) 


O 2 (^,0,0) 


= e m „ P ,e^n^ mn "0 Jpsfc 0^, 


(10.15) 


O 3 (0,0, 1p) 


^mnpq ^rstu 0^ ^ , 


(10.16) 




_ P P imnri ijpsk jlqtu 
c mnpq c rstu Y Y Y i 


(10.17) 




^mnpq ^rstu 'y' Vp^ ^ , 


(10.18) 




^mnpq ^rstu 1 0'"' ijj 


(10.19) 



An important observation is that all T's are vanishing for ( |10.8|) . Indeed, let us denote 
P j := UiF ij , then 

_ (10.20) 



^ijki = _ 2 ^\g\P [i u j] e kl . 
Then we straightforwardly find, for example, 

= 2|</| e mnw e rstM x wsfc P m u n {P l v? - P* u l ) e" e sk = 0. 



(10.21) 



This is zero because either the symmetric u n u q or symmetric P m pi is contracted with the 
antisymmetric e mnpq . Note that we on purpose write x~* psk as the second argument, because 



its form is arbitrary, not necessarily equal to ( [10.71 ). Analogously, we find: 

T 2 (iP, tP, X ) = A\g\ e mnpq e rstu P^ m vfte* P y uV fc X lqtu 

= 2\g\ e mnpq e rstu X lqtu P m u n (P^ vP - P p v?) e" e sk = 0. 
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10.22) 



It is a little bit more nontrivial to prove that T3 also vanishes. We have, explicitly: 

T 3 ( X , V, ■*!>) = M9\ £mn Pq erstu P [] u^e sk P9 u^e tu X mnn 

= \g\ e mnpq e rstu x mnn (P J u p P l u* - P p u>P l u q 

— P j u p P q u l + P p u j P q u l ) e sk e tu = 0. (10.23) 

The first and the last terms in the parentheses contain the symmetric combinations u p u q 
and P p P q which are vanishing when contracted with the antisymmetric e mnpq . The two 
remaining terms in the parentheses are reduced, by means of a relabeling of indices, to 
— P v vfi{P l vj — Pi u l ). Recalling that at the end we impose the symmetrization over the 
free indices k, I), we thus prove that T 3 (x, ip, ip) = 0. 

Since in all the three formulas ( |10.21| )-( |10.23D , the argument x^ psk is completely arbitrary, 
we can put x jpsk = ^ jpsk , in particular. Then ( |T0~2l| )-( [T0~2"^ ) yields that G(ip) = 0. As the 
next choice, we put x ipsk = ¥ psk ■ Then (|T0~T3l) combined with ( |IMIl )-( |RT2l ), yields 

Q{X) = Q{4>) + I (Ox + 2 + 3 ) . (10.24) 

It thus remains to compute the O-terms. The corresponding calculation is straightforward 
and simple if we use the representation ( |10.12j ). Then we find: 



Ox(0, ^, 0) = 2 (V>, 0, 0) = Os(0, 0, V) = - § sign(<?) - ¥ ij V fc|mn| ° L (10-25) 

Note that this result is valid for all possible tensors ip, not only (|10.8| ) which means that we 
can further use ( |10.25| ) for the future calculations involving more general nonlinear pieces 



(in particular, for the case with a nontrivial m lJ ). Using then ( |10.25| ) in ( |10.24| ), we get 



gv k \<j) + = — sign(s) yj\g\ g® 9 kl \ (10.26) 



_£ . / s r, r ?,..„. 1 
// 

Here we denoted 

gij ;= °gij _ V ^(i\mn\j) °g mn _ ( 1Q 27 ) 

y\d\ 

Hence, a purely dielectric nonlinear moving medium will in general exhibit the birefringence 
effect: the light will propagate in such a medium along the cone of the original optical metric 
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9 %] (one may call it ordinary ray), and along the second cone determined by the metric 9 lJ 
( "extraordinary" ray) . 

The explicit form of the extraordinary optical metric is obtained when we substitute 

( pmzq ) into ( [HT27D : 



9 l] + n e k{l u j) F kl u l - - u k e k(l F^ l u h 

e 



(10.28) 
(10.29) 



Special cases of this general result were considered in [[[J 

As an example, let us consider a medium in its rest frame. We use adapted coordinates 
such that u 
by 



8q = (l, 0). We additionally assume that the dielectric permittivity is given 



e = e + aE 2 . 



10.30) 



Here e and a are constant parameters. The components of the electric vector are defined as 
usual, E = E a = — Fo a . Then we find 



e 0a = _ e a0 = _ 2a E a^ ^ pj = ^ pij = E b\ 



(10.31) 



The spatial indices are lowered and raised with the help of the 3-metric S a b (we neglect the 
gravitational). As a result, from (|10.31|) we obtain 



P k e ki = (- 2a E 2 , o) , u k e ki = (o, - 2a E h ) . 



(10.32) 



Finally, using all this in ( |10.28| ) we find the components of the extraordinary optical metric: 



2a 



g ™ = ( 1 + E< 



^ab tab 



2a 



S ab E a E . 



(10.33) 



Here n := y/ejl is the refraction index of the medium. In this way we obtain a natural 
description of the optical Kerr effect (see [^,^8j, for example) when birefringence is induced 
by the applied electric field. 



20 



11. DISCUSSION AND CONCLUSION 



In this paper, we have performed a systematic analysis of the light propagation in the 
nonlinear electrodynamics on the basis of the Fresnel approach. We have considered two 
classes of models: (a) general nonlinear Lagrangian theories, and (b) moving nonlinear 
matter. In the former case, the Lagrangian ( |3.2| ) is an arbitrary function of the two electro- 
magnetic invariants, whereas in the latter case, the permeability and permittivity functions 
of the medium ( |10.1 ) depend arbitrarily on the electromagnetic field. 



The study of the first class of models reveals the generic nature of the birefringence effect: 
The quartic Fresnel surface reduces to the two light cones for all nonlinear Lagrangians. We 
show that the resulting optical metrics are always real and have the correct Lorentzian 



signature. In this way, we confirm and extend the recent results of fL3|J14 |. Furthermore, 
we are able to demonstrate the validity of the so called closure-no birefringence conjecture 
in the context of nonlinear electrodynamics: Birefringence is absent (and thus the quartic 
Fresnel surface reduces to a unique light cone) if and only if the effective constitutive (or 
jump) tensor satisfies closure property. 

The nonlinear moving matter with the constitutive law ( |10.1|) gives a sound example of 
a model in which the effective constitutive tensor naturally has nontrivial axion and skewon 
contributions. Accordingly, one should then expect that the Fresnel surface remains quartic, 



in general p3fl . However, for nonmagnetic material media, we show that birefringence is 
again the generic effect. The Fresnel surface factorizes into two light cones, one of which 
corresponds to the Gordon optical metric (independent of nonlinearities), whereas the other 
( |10.28| ) manifests the nonlinear properties of the model. The optical Kerr effect represents 
a particular example of our general derivations. 
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